Numerical study of the localization length critical index in a network model of 
plateau-plateau transitions in the quantum Hall effect 
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We calculate numerically the localization length critical index within the Chalker-Coddington 
(CC) model for plateau-plateau transitions in the quantum Hall effect. Lyapunov exponents have 
been calculated with relative errors on the order 10~ 3 . Such high precision was obtained by con- 
sidering the distribution of Lyapunov exponents for large ensembles of relatively short chains and 
calculating the ensemble average values. We analyze thoroughly finite size effects and find the 
localization length critical index v = 2.5f7 ± O.Of 8. 
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Plateau-plateau transitions in the quantum Hall effect 
have been one of the most challenging problems in con- 
densed matter physics during the last two decades. It 
is a interesting example of the localization-delocalization 
transition in two-dimensional disordered systems where 
a quantum critical point appears due to the breaking of 
time reversal symmetry. The most important problem 
in this area of research is the formulation of a quantum 
field theory describing the transition. The first sugges- 
tion in this respect appeared in Ref. [l], where the authors 
noticed that the presence of the topological term in the 
nonlinear sigma model formulation of the problem can 
result in the occurrence of delocalized states in strong 
magnetic fields. 

Later, Chalker and Coddington^ formulated a phe- 
nomenological model of quantum percolation based on 
a transfer matrix approach (refered to as the CC model 
hereafter). The numerical value 2.5 ± 0.5 of the crit- 
ical index of the Lyapunov exponent calculated within 
the CC model (see Ref. for a review) was in good 
agreement with the experimentally measured correlation 
length index v — 2.4 in the quantum Hall effect^ This 
success motivated considerable interest in the CC model 
and stimulated its further investigation until present 
days^i&£i&iiii2i22i In early studies 5 -^ the continuum 
limit of the CC model was related to supersymmetric 
spin chains, which was further developed by Zirnbauer, 9 
In Refs. [ll| and[l2| the continuum limit was also related 
to the conformal field theory of Wess-Zumino-Witten- 
Novikov (WZWN) type. Analyzing the representations 
of the PSL(2|2) conformal field theory they found one 
which gives a reasonable value of 16/7 ~ 2.286 for the 
correlation length index. Moreover, multi-critical scal- 
ing indices of the CC model were predicted to depend 
quadratically on the level q of multi-criticality within the 
WZWN model approach. 

Most intriguing developments in the plateau-plateau 
transition problem were reported later ; 13 ' 14 where the 
multi-critical behavior of the CC model was investi- 
gated. In both papers quartic dependence of the multi- 
critical indices of the parameter q was observed, in con- 
trast to earlier predictions ^ 1 ' 12 The latter suggested that 



the validity of the supersymmetric WZWN approach to 
plateau-plateau transitions in the quantum Hall effect 
is questionable. On the other hand, since the plateau- 
plateau transitions are of the second order, there is a 
conformal symmetry at the transition point and there 
should exist a conformal field theory describing it. Can- 
didate theories could be tested against the experimental 
data by comparing critical indices. Unfortunately, the 
precision of the available experimental indices is too low 
and does not enable us to identify the correct theory. 
Therefore, comparison to numerically calculated values 
can be more feasible for this task. However, reliable cal- 
culation of the localization length critical index with good 
precision has been known to be a very challenging task 
and there is still little consensus on the obtained values 
and especially their error bars. This paper is largely mo- 
tivated by the demand of such an accurate calculation. 

We have carried out numerical calculations of the 
Lyapunov exponent and the corresponding critical in- 
dex in the CC model taking finite size effects into ac- 
count. Our calculation strategy can be summarized as 
follows. Instead of finding the eigenvalues of the product 
T = Ylj=i Tj of a large number of transfer matrices Tj, 
we form a large ensemble of L r shorter products and ob- 
tain the distribution function of the smallest Lyapunov 
exponent 7 = (l/2L)lnA, where A is the eigenvalue of 
the matrix S = T^T closest to the unity from above. 
Then we calculate the mean value and the standard de- 
viation of the distribution. According to the central limit 
theorem the distribution of the mean value is normal 15 . 
This procedure allows us to effectively increase the prod- 
uct length up to L c ff — L r x L and reduce the error, 
which is on the order of 1/7/ £ e ff- The error can there- 
fore be reduced by collecting the statistics of eigenvalues. 
In finding the correct value of the critical exponent we 
take into account the finite size of matrices T and do the 
standard finite size scaling (FSS) analysis according to 
Refs. [llG]]GJ. The result we found for the correlation 
length index is v = 2.517 ± 0.018. 

The transfer matrix of the CC model consists of local 
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FIG. 1: Schematic view of the CC model network structure. 



2x2 matrices of the form 



WlR 

w 2 r 



l/t r/t 
r/t l/t 











W 2 L 



(1) 



where t = l/Vl + e 2x is the tunneling amplitude of elec- 
tronic states and r = vl— ir is the rotation at each node 
of the regular lattice shown in Fig. [T] The model param- 
eter x corresponds to the Fermi energy while phases on 
(i = 1, 2) are random variables from to 27r, arising from 
the randomness of the potential landscape. N of such 
local transfer matrices correspond to a vertical strip in 
Fig. Q] that forms 2N x 2 N transfer matrices W\ U\ and 
W2U2 for neighboring columns. The parameterization of 
transfer matrices W\ and W2 in neighboring columns is 
such that it is invariant under the rotation by 7r/2, the 
symmetry which is obvious from Fig. [TJ We therefore 
have two parameters x and x' corresponding to neigh- 
boring columns related as follows 



t = r' = 



1 



VT 



=, r = t' = 



1 



vT 



(2) 



The transition point x c is equal to 0. As the result, 
matrix elements of W\ and W2 are defined as follows: 

^l]2n+12n+l = Wl] 2n2n = V* > 

[Wi] 2n+12n = [W 1 ] 2n2n+1 = r/t, 
[W 2 ] 2n . 12n ^ = [W 2 ] 2n2n = l/t' , 
[W 2 ] 2n . 12n = [W 2 ] 2n2n _ 1 = r'/t' , 
[W 2 ] 1N = [W 2 ] Nl = r'/t' , 

where n = 1 . . . N and periodic boundary conditions are 
imposed on W2 ■ Matrices U have a simple diagonal form 
[U\,2\ nm — exp (ia n ) S nm . The final transfer matrix T 
can then be represented as a product Ylj=i W\U\jW 2 U 2 j. 

The largest matrix element of the transfer matrix T 
is of the order of r\ ~ (1 + e 2x ) L , growing exponentially 



with the system length L^L Because we are looking for 
eigenvalues of S = TT^ which are closest to unity, we 
need to compute S with the precision of about 2 log 10 t?, 
otherwise numbers of the order of unity are unreliable 
due to round-off errors. We used arbitrary precision cal- 
culations with adaptive precison increasing it gradually 
for each matrix product. We calculated the transfer ma- 
trix T of the system of length L on a grid {N, x} using 
a graph-like algorithm. Instead of computing product of 
local transfer matrices sequentially, we calculate it on a 
full fc-way tree with adaptive precision. At each node we 
estimated the required precision and set the involved ma- 
trices to it before multiplying them. In all calculation we 
used k = 7, which allowed us to perform the vast majority 
of the matrix multiplications (at the leaf node level) with 
the machine precision. When the matrix S = TT^ was 
computed we used the shift-and-invert Lanczos algorithm 
to calculate a set of selected eigenvalues of S which are 
closest to unity. Typically, we calculated six eigenvalues 
with the accuracy of 10~ 5 , checking that these eigenval- 
ues appeared in pairs and satisfied log \2i — — log X21+1 
(where i = 0, 1, 2 and Ai is the eigenvalue closest to unity 
from above, i. e. Ai > 1). We then calculated Lyapunov 
exponents — (1/2N) log[Aj(iV, a;)]. Hereafter, we only fo- 
cus on the smallest Lyapunov exponent ^n^)- 

The tree-like algorithm provides access to partial trans- 
fer matrices and allows us to calculate Lyapunov expo- 
nents at different node levels, not only at the root node. 
The latter offers the possibility to collect and study the 
statistics of Lyapunov exponents. We analyzed distri- 
butions of 7Ar(a;) for different system lengths L. Fig- 
ure[2ja) displays the mean value 7jv(^) of the distribution 
and its estimated standard errors calculated for N = 60, 
x = 0.404, and L = 343, 2401, 16807. Larger error 
bars show the standard deviation <jn{x) of the distribu- 
tion at each node (i.e., at each L). Smaller error bars 
represent the standard error of the mean value itself, es- 
timated as 



(3) 



a N (x) = a N (x)/\/L 



where L r is the size of the Lyapunov exponent ensemble 
at each node. Figure [2](b) demonstrates that the stan- 
dard deviation <Jjsi(x) varies as 1/y/L with the system 
length, confirming the applicability of the central limit 
theorem and, hence, the definition ([3]). 

Figure [5] shows that values of the disorder-averaged 
mean Lyapunov exponent for N > 2401 are the same 
within their error bars. We calculated these quantities 
for different sets of N and well as for longer chains, 
and found the same result. Therefore, rather than calcu- 
lating the exponent for a very long chain, we can calculate 
the mean Lyapunov exponent for a large number L r of 
relatively short chains, which is more efficient. Varying 
L r we can achieve any given target accuracy; in our cal- 
culation we used L = 4802 with the typical ensemble size 
of N r w 2000. 

Having calculated disorder-averaged Lyapunov expo- 
nents, we use the standard FSS analysis of the data, for- 
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FIG. 2: a) Mean values of the 7jv (^-distribution and esti- 
mated error bars calculated for N — 60, x — 0.404 at dif- 
ferent node levels (L = 343, 2401, . . . , 16807). The horizontal 
dashed line is a guide for the eye, showing the mean value for 
L = 16807. Larger error bars show the standard deviation 
of the distributions (tn(x), while smaller error bars are equal 
to gm(x) and estimate the error of the mean value, b) Full 
circles show the standard deviation <jn{x) as a function of the 
system length L and the solid line is the best 1 / \[~L fit to the 
data. 



mulated in Ref. 16 and extended in the Refs. ll7lll8lll9l . in 
order to obtain the localization length index v. The Lya- 
punov exponent is believed to have scaling behavior (see 
for example Refs. [l?] andl20l and references therein), and 
finite size effects can be accounted for by the following 
formula for the scaling function T(N,x) which approxi- 
mates Njn(x) in the vicinity of the critical point 



T(N,x) = F {N 1/,J u ) + Fi{N 1/v u Q )N y u x 



(4) 



where Fq(-), Fx(-), uq(-) and iti(-) are universal, indepen- 
dent of L functions. The first function Fq(-) is the con- 
tribution of the main operator in the corresponding con- 
formal field theory, which defines the correlation length. 
The second function -Fi(-) results from the operator with 
the anomalous dimension which is close to that of the 
main one. We used the following formula to fit to the 
data« 



T(N,x) 



3 



, AO 2r. 



n=0 



1 

u^N 1 '" +a loUl {x)Ny (5) 



where u (x) = x + b 3 x 3 and Ui(x) — 1 + c 2 x 2 , so the 
series for F Q , F\, u Q , and Ui were truncated at orders 6, 
0, 3 and 2, respectively. Further increase of these orders 
turned up to be impractical as it deteriorated the result. 

Because all mean Lyapunov exponents have different 
error bars we used a weighted fit with [iV(7L(:r)]~ 2 as 
weights. We also put the following constraints on the 
parameters: 0.5 < a 00 < 1, |fflo2n|, Ko|, |c 2 | < 2, 
2.1 < v < 2.7 and -0.5 > y > 0. We checked that 
the result is not affected by the constraints while they 
stabilize the procedure and improve the convergence of 
the nonlinear fit. All fits which converged to a constraint 
boundary were discarded. Global optimization methods 
such as the simulated annealing were used to find the 
best fit. The result of such a fit is presented in Fig. [3] 
All curves intersect at x rs 0.025; note that the ordering 
of the curves at x — suggests an additional constraint, 
namely a\Q > 0. 

In order to estimate error bars of the parameters we 
used the standard resampling technique ~i Up to 1000 
generated synthetic data sets of mean Lyapunov expo- 
nents were drawn from the corresponding normal dis- 
tributions centered at 7at(x) and having the standard 
deviation aN{x). Fitting to these data we obtained dis- 
tributions of the parameters and calculated their mean 
values and standard errors. The following parameter 
values were found: for the case c\ / we obtained 
v = 2.519 ±0.023, y = —0.32 and the typical goodness of 
fit parameter was \ 2 — 864.5 for 107 points. For c\ = 0, 
v = 2.517 ± 0.018, y w -0.24 and X 2 = 848.2. We note 
that the parameter y is always small, which suggests that 
finite size effects are very pronounced in the CC model. 

The curves for Lyapunov exponents for different sizes 
N are intersecting in a region close to x = 0. The crossing 
point of a pair of such curves corresponding to consecu- 
tive values of N shifts towards the origin on increasing N. 
The latter is the consequence of the condition <2iq > 0, 
which increases the role of the next to the leading op- 
erator in the problem for small N. This fact stresses 
the importance of the second operator together with the 
main one in the analysis of possible candidate conformal 
field theories which describe plateau-plateau transitions 
in the quantum Hall effect. 

The obtained value of the critical index v = 2.517 ± 
0.018 disagrees with early results^ ' 22 i 23 which could 
be attributed to the quality of the data, system sizes 
reached, and finite size effects. Those effects prooved 
to be extremely important for accurate calculation of the 
critical exponent even in the standard Anderson modelj^ 
in which case the irrelevant exponent y ~ —3 is much 
larger than in the CC model (finite size effects are there- 
fore less pronounced) . On the other hand, our critical ex- 
ponent is smaller than the one calculated by Slevin and 
Ohtsuki, 20 who obtained v = 2.593 ± 0.006. However, 
both high precision numerical results are considerably 
above the experimental value of v = 2.38±0.06 measured 
recently in GaAs-AlGaAs heterostructures i 24 i 25 The lat- 
ter fact emphasizes the necessity to investigate the valid- 
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FIG. 3: Simpulated data and the fit of Eq.[5]to them (ci = 0). 
From bottom to top (for x > 0.2): N = 10, 20, ... , 60. 

ity of the CC model further. 

In summary, the obtained critical exponent v > 2.5 
suggests that the rational value 7/3 which is in agree- 
ment with some early calculations, could be questioned. 



Our technique is based on the calculation of the prod- 
uct of transfer matrices preserving the required precision 
during iterations and collecting large statistics to achieve 
a given target precision. Finite-size effects turn out to be 
very pronounced within the framework of the CC model 
and should be taken into account properly in order to ob- 
tain reliable results. In particular, the calculated value 
of the critical exponent suggests that some WZWN-type 
models based on the conformal field theory should be 
reconsidered, which demands new developments and ap- 
proaches in the formulation of the continuum limit of the 
CC model as well as further studies and more accurate 
calculations of Lyapunov indices. 
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